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Abstract 

The n + 1-dimensional Milne Universe with extra free directions is used to construct 
simple FRW cosmological string models in four dimensions, describing expansion in the 
presence of matter with p = up, k = (4 — n)/3n. We then consider the n = 2 case 
and make SL(2, Z) orbifold identifications. The model is surprisingly related to the null 
orbifold with an extra reflection generator. The study of the string spectrum involves the 
theory of harmonic functions in the fundamental domain of SL(2, Z). In particular, from 
this theory one can deduce a bound for the energy gap and the fact that there are an 
infinite number of excitations with a finite degeneracy. We discuss the structure of wave 
functions and give examples of physical winding states becoming light near the singularity. 
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1. Introduction 



Deep questions in cosmology, such as the nature of the big bang singularity or initial 
boundary conditions, need to be addressed in the context of a quantum gravity theory. 
Understanding string theory in time-dependent backgrounds leads to a number of technical 
problems which have been considered in many recent works [|T1-[T6||. In particular, in || 



the null orbifold model ||17|1 - a time dependent geometry consisting in flat space with 
some identifications - was studied in detail. The string model of this Lorentzian orbifold 
contains some instabilities which are reflected in IR divergences in scattering amplitudes. 
Their physical origin is a large back reaction in the geometry as soon as particles are 
introduced 



In the present work we will be interested in the general class of time-dependent locally 
flat spacetimes obtained from the n + 1 dimensional Milne Universe. They can be viewed 
as a particular case of hyperbolic compact ificat ion in string theory (see ||18|| ), which have 
recently attracted some interest as they lead to interesting cosmologies |IT^] . We will 
consider in more detail the n = 2 case, which, with proper identifications, is related to the 
null orbifold. Here we make orbifold identifications on the upper half plane by using the 
(elliptic) modular group SL(2, Z). 



Milne spaces in the context of inflationary cosmology were studied in [20]. String 



models in 1+1 dimensional Milne space were discussed in [|T7| ,p],|3|fl . Discussions on higher 



dimensional Milne spaces can be found in fll^Jlql , and more recently in [|16 



In section 2 we discuss the simplest string cosmology based on flat Milne space. We will 
see that already these simple backgrounds can be used to obtain interesting cosmologies 
in four dimensions. This section makes a connection with recent work on hyperbolic 
compactification ||19|| , and it can be viewed as a physical motivation for string models 
based on Milne space. In section 3 we introduce the SX(2, Z) orbifold model. The study of 
this model involves harmonic analysis on symmetric spaces and number theory, from which 
we will extract some results for the string spectrum. We also discuss the construction of 
wave functions and physical string states. 

2. String cosmology from n + 1-dimensional Milne Universe 

The n + 1-dimensional Milne space is described by the metric 

ds 2 = -dt 2 + t 2 dHl , (2.1) 
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where dH 2 is the arc element of the hyperboloid or upper half n plane. In Poincare 
coordinates, it is 

dH* = ^{dz 2 + d4 + ... + dzl_ 1 ) . (2.2) 
The space is flat, as it is evident upon introducing Cartesian coordinates as follows: 

U=-, V = tz + -Y ( zt, X l = ^. (2.3) 
z z -f— : ' z 

This provides the embedding of the hyperboloid in n + 1 Minkowski space, 

ds 2 = -dUdV + dXf , 

where the hyperboloid is described by 

t 2 = UV-X 2 -...-X 2 _ 1 , (2.4) 

which exhibits the SO(l,n) isometry of H n . 

Alternatively, one can use angular or Milne coordinates 

dH 2 = dp 2 + sinh 2 p dfl 2 ^ ; (2.5) 

and define 

T = icoshp, i? = isinhp, (2.6) 
in terms of which the metric ( [2.1| ) reads 

ds 2 = -dT 2 + dR 2 + R 2 dSl 2 n _ x . (2.7) 

Now let us consider type II string theory in the following background 

ds\ Q = -dt 2 + t 2 dHl + dx\ + dx\ + dx\ + dy\ + ... + dy 2 _ n , (2.8) 

where the y coordinates describe compact internal dimensions. All other supergravity 
fields are trivial. This is an exact solution (to all alpha prime orders) of string theory, 
since the Riemann tensor identically vanishes (though it may receive quantum string-loop 
corrections). Similarly, one can write down the analogous eleven dimensional metric, which 
is a solution of M theory. The internal space described by the y coordinates can be replaced 
by any Ricci flat space, giving a more general class of cosmological backgrounds which are 
solutions to leading order in a' . 
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It is surprising that an interesting four dimensional FRW cosmology can be obtained 
from the model (|2.8|). First, we replace the hyperboloid H n by a finite volume space H n /T, 
where T is a discrete subgroup of SO(l,n) such that the space has finite volume, as in 
T^,|19[. Then we compactify to four dimensions. To obtain the four dimensional Einstein 



frame metric, we write the ten-dimensional metric in the form 

ds 2 = e 2a ^ds 2 E + e 2b ^dH 2 + dy 2 + ... + dy 2 _ n , (2.9) 

ds\ E = e 2c{t) (-dt 2 + dx\ + dx\ + dx\) , 

and the condition that ds 2 E is in the Einstein frame is e 2a e nb = 1. Comparing to ( |2.8|) , 
one obtains 

ds\ E = t n (-dt 2 + dx\ + dx\ + dx\) , (2.10) 

or 

ds\ E = -dr 2 + r^{dx\ + dx\ + dx\) . (2.11) 

This corresponds to 4d Einstein equations coupled to an energy momentum tensor of a 
perfect fluid with 

4 — n 

P = K P , K= — — . (2.12) 

in 

Although we have started with vacuum Einstein equations in ten dimensions, the four 
dimensional Einstein metric describes a homogeneous and isotropic space in presence of 
matter. This matter is of course the scalar field associated with the modulus representing 
the volume of the hyperbolic space. Interestingly, the above metric is the asymptotic (large 
time) form of the models of [|I9" |. For n = 4, it describes a universe filled with dust, and 



for n > 4 a universe filled with negative pressure matter. 

Since the models are based on a flat ten dimensional geometry, the n + 1 dimensional 
Milne Universes provide a simple setup for the study of interesting cosmological string 
models. 

So let us consider strings propagating in this space. An important question is whether 
the string model is exactly solvable. To start with, consider the model ( |2.8| ) based on H n 
with no identification, i.e. T = 1. The string equations of motion in Cartesian coordinates 
( |2.3|) are solved explicitly in terms of right and left moving modes. However, from the 
relation ( |2.4|) it follows that 

UV-X 2 >0. (2.13) 
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If the physical space is restricted only to this Milne patch, say with t > 0, then the 
string coordinates are subject to the constraint that the string lives in the interior of the 
future directed light cone; the space is not geodesically complete and a full description 
requires boundary conditions at the light cone surface. In string theory, it is possible that 
consistency also requires the inclusion of the past light cone. In this case, the geometry 
would describe a universe contracting to a big crunch which makes a transition to an 
expanding big bang universe. 

It is non-trivial to impose the condition ( |2.13|) in string theory. On the other hand, if 
the full space U, V, Xi is considered, closed timelike curves can arise in the exterior of light 
cone as a result of identifications. 

We briefly comment on a possible approach with a boundary on the light cone surface. 
A string which at a given instant fully lies in the interior of the future light cone, will satisfy 
the condition UV — Xf > at all future times. Causal propagation implies that all closed 
strings in the physical space emerge from the light cone surface, so there is some interval 
of time where the condition UV — Xf > is satisfied only for a piece of the string (an 
example is given at the end of section 3). It is natural to describe the strings intersecting 
the light cone in terms of open strings with ends attached to the null light-cone surface, 
which can be thought of as a "null D brane" . In this picture, the set of open string states 
attached to the null surface determine the initial state, i.e. they would be described by 
some "in" state in the open string Hilbert space. The closed strings living in the interior of 
the light cone would originate by emission from the null D brane. Thus the time evolution 
of the cosmology would be determined by the initial state on the null D brane and the 
decay properties of this initial state. 

Let us now comment on the supersymmetry of these models. The hyperbolic upper 
half n-plane H n has Killing spinors transforming in the spinorial representation of SO(l, n— 
1). They were explicitly constructed in [^TJ (see also discussions in |2^ , |T8| , [23[1 ) . The simplest 



model (|2.8| ) with V = 1 is supersymmetric. When T is not trivial and the space H n /T has 
finite volume, it was shown in |H| that for n =even supersymmetries are always broken 
by the identifications. For n =odd and n > 5, the analysis of does not exclude that an 



appropriate choice of V could give a supersymmetric model with finite volume hyperbolic 
space. It would be interesting to find some example. 



3. SL(2, Z) orbifold 



Let us consider the case n = 2 and recall its relation to the null orbifold ||17|1 . The 
metric is 

t 2 

ds 2 = -dt 2 + — (dx 2 + dy 2 ) . (3.1) 

The constant time slices describe the hyperbolic upper half plane H 2 . Introducing U = t/y, 
V = ty, we obtain 

ds 2 = -dUdV + U 2 dx 2 . (3.2) 

This form exhibits an orbifold singularity at U = moving at the speed of light. Intro- 
ducing new coordinates 

u = U , v = V + Ux 2 , X = Ux , (3.3) 

we get 

ds 2 = -dudv + dX 2 , u = X -X 1 , v = X + X 1 . (3.4) 



The null (or "parabolic" ) orbifold geometry of [TJ] , studied by Liu, Moore and Seiberg , 
is obtained by orbifolding the isometry x = x + 1. In the coordinates (m, v, X) this leads 
to the identification 

(u,v,X) = (u,v + un 2 + 2Xn,X + n u) (3.5) 

Now we consider another geometry: we replace H2 in ( |3.1| ) by T = SL{2, Z) \ H2M 
SL(2, Z) transformations are generated by z —>■ z + 1 and z — > —1/z. The first transfor- 
mation is precisely what leads to the null orbifold (|3.5|) . Now we have the extra orbifold 
identification z — > — 1/z, with z = x + iy (and t — > t). This isometry is surprisingly simple 
in the coordinates (m, v,X). It just corresponds to the spatial reflection: 

u — > v , v — > u , X — > —X , (3.6) 

where we have used the relation to the original coordinates (t, x, y) given by 

u = -, v= t -(x 2 + y 2 ), X=-. (3.7) 

y y y 



1 In the notation of Q one has x + = u, x = v/2. 

2 A brief discussion of this geometry is in section 5 of [24|. 



5 



It should be noted that it is not sufficient to add this orbifold identification ( |3.6|) to 
the string model of || in order to have the fundamental domain T as constant time 
slices. In addition, one has the restriction uv — X 2 > 0. So one can distinguish two 
related but different models: the "null orbifold with reflection" or just U SL(2, Z) orbifold", 
obtained simply by adding the reflection identification (|3.6| ) to the null orbifold model of 
and the "restricted SX(2,Z) orbifold model" (|3.1|) having the finite volume space 



T as constant time slices. Both models arise as a quotient by the <SX(2, Z) subgroup 
of the Poincare isometry group 5'L(2,R) and have SL(2, Z) symmetry. The former has 
another patch uv — X 2 < (i.e. UV < 0), where the metric (|3.2| ) can be put into 



2 2 

the form ds 2 = dt 2 + ^{dx 2 — dy 2 ) = dt 2 + ^dz + dz_ , with U = —t/y, V = ty and 



z± = x ±y. In this patch, the identification ( |3.6| ) corresponds to the isometry z + — > l/z+, 
Z- — > 1/z-, t — > t. In the null orbifold model with reflection it seems more convenient to 
work in (u, v, X) coordinates. This model has closed time-like curves arising in the region 
uv — X 2 < 0. To see this, one can consider the image of the point u = v = 0, X = 1 
under an SL(2, Z) transformation with parameters {a, 6, c, d}, ad— be = 1.1 For be > 0, the 
images («', v', X') = (2ab, 2cd, l + 26c) are time-like separated from the original point. This 
transformation is a combination of reflection (|3.6| ) and the parabolic orbifold identification 
( |3.5|) . The identification ( |3.6| ) alone, or the identification ( |3.5| ) alone, do not produce CTC. 

Here we shall use properties of modular functions in T so in the remainder we will 
consider the restricted SL(2, Z) orbifold model. The two string models should be closely 
related. 

In addition to the singularity at U = t/y = 0, the geometry contains the singularities 
of T at \z\ = 1, x = ±1/2 for all t. In the flat coordinates, it corresponds to the subspace 

u = v = ±2X . 



There are no surviving supersymmetries, they are broken by the SL(2, Z) identifications. 

Let us now investigate the string spectrum of this model and the structure of vertex 
operators. To leading order in a', physical modes of the string spectrum obey, in a suitable 
gauge, the Klein-Gordon equation. This applies to massless and massive scalar fluctuations, 
as well as to non-zero helicity components of a massive or massless mode. The only 
exception are winding modes, which cannot be described in terms of local fields of the low 

3 We thank Hong Liu for providing this example. 
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energy string effective action. We will construct examples of invariant winding states at 
the end. 

First, consider a massless scalar field propagating in this geometry. To exhibit the 
SL{2, Z) invariance, it is convenient to work in (£, x, y) coordinate system. We set the 
transverse momentum components pi = 0. The Laplace equation A3 1 ]/ = is 

A 2 tf = d t (t 2 d t V) , (3.8) 

2 ^ dx 2 dy 2 
We can solve it by separation of variables, \1/ = i/j s (x, y)(p s (t), with 

A 2 V> S = s(s - 1)V> S • (3.9) 

This implies that the solution is given in terms of Maass waveforms. These are modular 
functions which are eigenfunctions of the Laplace operator A 2 with at most polynomial 



growth at y — > 00 [25 



Let us first obtain the wave function by simply starting with the flat solution and 
summing over images. In the coordinates (u, v, X) the solutions of the wave equation 
A3 1 !/ = are plane waves with momentum {p u ,Pv,Px), with 4p u p v = p 2 x . Then a invariant 
wave function would have the generic form 



* = X) / dpudp v dp x <p(p u ,p v ,p x ) e*»«+ < P»»+ i - 



ip x X 



images 



J dp u dp v (p(p Ul p v ) ^2 ' e l y pu \( a+cx)z+b+dx \ , X = y/p v /Pu , 



(3.10) 



a,b,c,d 



where Yl' stands for a sum over integers a,b,c,d obeying ad — be = 1, and we have 
used the relation ( |3.7| ). Normalizability of the wave function ( |3.10| ) (with the invariant 
measure dxdy/y 2 ) is a delicate issue. To illustrate this point, consider plane waves with 
[PuiPv] = q(m 2 ,n 2 ), with px again fixed by the mass shell condition 4p u p v = p 2 x . They 
are of the form 

®=[dq c(q) e iq i lm+nzl2 . (3.11) 

(m,n)9*(0,0) 

Now consider a c(q) of the form: c(q) = c s q s ~ 1 9(q), with s being a complex number and 
9(q) is the step function. The integral in q is a Mellin transform. Computing this integral, 
one obtains 

* = 2c s T(s)C(2s) t~ s E s {z) , (3.12) 
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where E s is the non-holomorphic Eisenstein series 



E s = Y i —rr ■ 3.13 

2C 2s ^ m + nz 2s y ' 

SV ' (m,n)/(0,0) 1 1 



It obeys the functional relation E s = const i?i_ s , so with no loss of generality one can 
restrict to Re s > 1/2.0 At y = oo, E s behaves as y s . Therefore none of these wave 



functions (|3.12| ) are normalizable. 

Normalizable wave functions can be written in terms of cusp forms, which are auto- 
morphic functions that decrease exponentially at infinity. They are bounded in H2. Let 
us summarize some known facts, referring to |3S[ for details: 

i) If Re s > 1/2 and s 7^ 1 , then the vector space V s of Maass waveforms with eigenvalue 
s(s — 1) is generated by the non-holomorphic Eisenstein series E s . For s = 1, the vector 
space V s is generated by the constant. 

ii) If Re s = 1/2, the vector space V s is generated by E s and the cusp forms v n , n > 1. 

iii) The vector space of cusp forms v n is different from {0} only for eigenvalues s(s — 1) 
with Re s = 1/2, for an infinite number of values s. The set of values is discrete. However, 
an explicit basis is not known, nor an explicit analytic example of a cusp form (they can 
be easily found numerically). 

Thus cusp forms constitute the discrete part of the spectrum, while E s constitute a 
continuous part of the spectrum. 

A general eigenfunction in L 2 (SL(2, Z) \ H) is given in terms of the Roelcke-Selberg 
expansion, which is the analog of a Fourier expansion in T . Define the inner product in 
the standard way as (f,g) = J :F g(z)f(z)y~ 2 dxdy. Then any if; in L 2 (SX(2,Z) \ H) has 
the expansion 

ip{z) = J2(i>,v n )v n + -^ [ &,E s )E s (z)ds , (3.14) 

where vo = y3/V- 

Any cusp form is part of a vertex operator of the string model having non-zero quan- 
tum numbers in T. A similar discussion applies for the wave function of a massive scalar 
field: we solve the equation 

(A 3 -M 2 )* = 0, (3.15) 



4 E s has a pole at s = 1. E s with s = 1/2 has a zero which is cancelled in (3.12) by the factor 
C(l), leading to a y 1 ^ 2 logy behavior at y = 00. 
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by separation of variables, \& = i(j s (x, y)<p s {t). The function ifi(x, y) is then an eigenfunction 
of the Laplace operator on J 7 , so it is also given in terms of cusp forms. Using eq. (|3.9|), 
<p s is then determined by 

[d t t 2 d t - s(s - 1) + M 2 t 2 ]<f) s {t) = • (3.16) 

Hence <f> a (t) = J ±(s _ ^{Mt)/y/t. 

The cusp form can be expanded in terms of Bessel functions as follows: 

i/) s (z) = CLn y 1 ^ 2 K s _ l (27r|n|y) exp(27rmx) , s = — \-ib , 6gR, (3-17) 
— ' 2 2 



n/0 



where a n are constrained by eq. ( |3.9| ). Alternatively, one can make calculations in the 
covering space and compute S-matrix elements by using plane-wave vertex operators of flat 
space and adding the images. However, we have seen above that naive vertices constructed 
by a sum over images as in ( |3.iq ) may not be normalizable (e.g. c(q) = J2kLo Ck< ^ does 
not give a normalizable wave function ( |3.11| )). SL(2, Z) invariance severly constrains the 
form of vertices and amplitudes. In what follows we shall deduce a few interesting facts 
about the spectrum using theorems about cusp forms. 

The ten dimensional string theory based on the SL(2, Z) orbifold has additional free 
directions (xi, ...,xj). For a massless scalar fluctuation obeying Ai \I/ = 0, with no string 
oscillations in the t, x, y directions, the vertex is of the form 

$ = J d 7 p o(p) <f>,(t) ip a (x,y) dx.dx e ip ^ , (3.18) 



where <f> s (t) is determined by (|3.16| ) with M 2 — > M 2 + p 2 (in this case, M 2 = 0). In the 



frame Pi = 0, the equation is of the form dtt 2 dt<j) e — e<p e . The parameter e, being the 
eigenvalue of the time-derivative part of the Laplace operator, measures the energy scale 
of the Kaluza-Klein excitation (bearing in mind that in general energy is not conserved 
since the background is not static) .1 The spectrum formula is just given by 

e = s(s- 1) , (3.19) 

and 4> s (t) = a s t~ s + b s t s ~ 1 . For s = we have e = 0, corresponding to a state with trivial 
quantum numbers in J 7 , being the only normalizable state with Re s ^ 1/2. For other 



For massless excitations, the energy parameter e is associated with the scaling t — > At. 
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states, normalizability requires that s = 1/2 + ib, so that e = — 1/4 — b 2 is negative definite. 



Which is the value of e for the first non-trivial excitation? A theorem |25[ establishes that 
s(s — 1) < — 3n 2 /2. This implies that the energy e of the first Kaluza-Klein state obeys the 
bound e < — 37r 2 /2. Precise numerical estimates of the first and the next eigenvalues are 



summarized in |25[ . Because the spectrum contains both a discrete as well as a continuum 
part, one should not expect that at low energy the theory will look D — 2 dimensional. 

Another important theorem establishes that the vector space of cusp forms with a 
given eigenvalue is finite dimensional. This means that there is a finite number of normal- 
izable excitations at each energy level e. 

Finally, property iii) implies that the levels for cusp forms are discrete and that the 
total number of normalizable excitations is infinite. There is also an asymptotic formula 
from which one can estimate the total number iV(eo) of normalizable excitations above an 
energy level eo (or with \e\ < |eo|), counted with multiplicity. For large |eo|, one finds the 
remarkable formula N(eo) = |eo|/12. This is a consequence of the Selberg trace formula 

& 

Similar results apply to quantum numbers of massive string modes. 

Let us now construct physical winding states associated with the twisted sector of the 
orbifold identifications ( |3.5|) , ( |3.6|) . In the coordinates (u,v,X), the solution to the string 
equations are expressed in terms of free right and left moving oscillators. We only need 
to impose boundary conditions and solve the constraints. We first consider the winding 
solution (cr± = a ± r) 

u = p v r, X = ^(al -a 2 ) = ^ar , X(a + 2n) = X(a) + n u . (3.20) 
87T 2tt 

Then v = v + (a + ) + t>_(a_) is determined by solving the constraint equations T ++ = 
T__ = 0. This gives 



n Pv I 3 3\ n Pv I 3 i o 2 \ /, 01 \ 



2 2 
"^y_ i 3 _ 3 \ _ n 1 

24n 2 [a+ °- ] ~ 12tt 2 

,2„ 



v(a + 2tt) = v(a) + n 2 u(a) + 2nX(a) . 
In terms of (t, x, y) coordinates, the solution takes the form 
, n Pv 2 na n 

t = — =— r , x = — , y = — 7 =- r . (3.22) 
2V3n 2tt 2V3tv V ; 
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For this state to be SL(2, Z) invariant in H2, one adds all the images in the upper half 
plane. 

Next, consider the following solution which is a winding state for the twisted sector 
of©: 

X = 2L sin(a/2) cos(r/2) , X 1 = 2L cos(a/2) cos(r/2) , (3.23) 
u = X - Xl , v = X + X 1 . 
From the constraints, one gets Xq = Lt. This solution satisfies the boundary condition 

u(a + 2tt) = v(a) , v(a + 2n) = u(a) , X(a + 2n) = -X(a) . 

In the covering space, it describes a pulsating circular string loop, which contract to zero 
at r = 7T and has maximum radius at r = 0, 2n, etc. 

In the (£, x, y) coordinate system, the solution becomes 

t = Ly/r 2 -4 cos 2 (r/2) , 

= 2sin(q/2)cos(r/2) = y/g - 4 cos 2 (r/2) 

r- 2 cos(a/2) cos(r/2) ' y r - 2 cos(a/2) cos(r/2) ' 1 ' ' 

with the addition of the images. Note that in these coordinates the solution exists after a 
r = 1.48 such that r 2 - 4cos 2 (r/2) > 0. For r < r , the condition uv - X 2 > is not 
satisfied for all points of the string, i.e. part of the string is in the patch uv — X 2 < which 
is not covered by the (t, x, y) coordinates. Different points of the string pass through the 
singularity u = at different times X = Lt. As r gradually increases and goes over r , 
the string fully enters into the region uv — X 2 > 0. At this point, the string has already 
left the singularity u = 0. In the upper half plane, when r is slightly above tq, the string 
describes a long loop. When r = it, it contracts to a point at x = 0, y = 1. Then it keeps 
pulsating around this point, with a maximum radius that decreases with time. 

In the Hilbert space, these classical string configurations ( |3.22| ), (|3.24j ) are described 
by physical string states with large quantum numbers, so that the classical description 
applies. This requires large p v in the case of ( |3.20| ), and large L - which implies large 
occupation number of the lowest frequency string oscillator - in the case of ( |3.23j ), ( |3.24|) . 
The state (|3.22|) which winds around x becomes "massless" near U = (the light-cone 
energy p u is proportional to U 2 n 2 /p v , see ( |3.2| ), (|3.3p). This is a reflection of the singularity 
of the geometry at U = 0. 
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The instabilities pointed out in [12] may be absent in the <SX(2, Z) orbifold model. 
Indeed, this instability originated from the gravitational interaction of plane waves and 
their images. Although here the spectrum has also a continuum part which may lead to 
wave interactions, this part is severely restricted by SL(2, Z) symmetry, so the argument 



of [jy! does not seem to directly apply to this model. The states in the discrete part of the 
spectrum have finite motion. The corresponding wave functions are regular (bounded) on 
J 7 , though they exhibit a singular behavior at t = 0. It would be interesting to investigate 
the partition function and higher point scattering amplitudes to establish if the string 
model is regular or singular. In the SL(2, Z) orbifold model, the extra reflection symmetry 



( plf ) adds new terms to the partition function computed in ||. 
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